Analytical Potential Energy Function for the 
Ground State XiS+ of LaCl 



Lin-Hong Chen*, Ren-Cheng Shang 

Center for Astrophysics, Tsinghua University, Beijing 100084, China 
Abstract 

The equilibrium geometry, harmonic frequency and dissociation energy of lanth- 
anum monochloride have been calculated at B3LYP, MP2, QCISD(T) levels with 
energy-consistent relativistic effective core potentials. The possible electronic state 
and reasonable dissociation limit for the ground state are determined based on 
atomic and molecular reaction statics. Potential energy curve scans for the ground 
state X^Ti~^ have been carried out with B3LYP and QCISD(T) methods due to 
their better performance in bond energy calculations. We find the potential energy 
calculated with QCISD(T) method is about 0.5 eV larger than dissociation energy 
when the diatomic distance is as large as 0.8 nm. The problem that single-reference 
ab initio methods don't meet dissociation limit during calculations of lanthanide 
heavy-metal elements is analyzed. We propose the calculation scheme to derive 
analytical Murrell-Sorbie potential energy function and Dunham expansion at equi- 
librium position. Spectroscopic constants got by standard Dunham treatment are 
in good agreement with results of rotational analyses on spectroscopic experiments. 
The analytical function is of much realistic importance since it is possible to be 
applied to predict fine transitional structure and study reaction dynamic process. 
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1 Introduction 



The transition metal monohydrides and monohalides have been extensively 
studied over several decades because they are of considerable interest in various 
fields such as astrophysics, catalytic chemistry, high-temperature chemistry 
and surface material [1-3]. In particular diatomic molecules of La are of great 
importance as test cases for modeling the role that d electrons play in chemical 
bonds, since lanthanum has the simplest open d shell electronic configuration 
[(core)5(i^6s^]. There have been significant efforts made in both spectroscopic 
investigations and ab initio calculations in order to get a better understanding 
of the electronic structure of such molecules. The visible band systems of 
monohalides LaH, LaF, LaCl are observed and rotationally analyzed [4-8]. A 
number of theoretical calculations also exist using ligand field theory and ab 
initio methods at all levels[9-12]. 

However up to now theoretical studies mainly focus on molecular equi- 
librium structure. Only several papers were devoted to calculations of global 
potential energy curves of LaH and LaF[13-15], which is not compatible with 
so many experimental spectroscopic data. In the case of LaCl, there is even 
no report on theoretical calculation of its structure and molecular properties. 
Global potential energy curve is significant for the investigation of collision 
reaction dynamics when these compounds are used as catalyst. In the present 
paper potential energy curve and analytical potential energy function, as well 
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as some molecular properties for LaCl molecule are firstly calculated and pre- 
sented with energy-consistent Relativistic Effective Core Potentials (RECPs) 
and the corresponding valence basis sets including diffuse and polarization 
functions. 

2 Calculation Method 

Quantum chemical calculations of lanthanide complexes are relatively dif- 
ficult due to complex electronic structure and large systems. A large number 
of low-lying electronic states of different spatial and spin symmetries may 
mix heavily. Electron correlation and relativistic effects can be very consid- 
erable, and spin-orbit interactions are often quite important for these sys- 
tems. Accurate theoretical calculations have to take account of all these ef- 
fects properly. As for the ground state X^T,'^ of LaCl, first-order spin-orbit 
interactions needn't to be considered. Although in recent years some advance 
has been achieved in four-component all-electron relativistic calculations [10- 
12], they are too costly. We adopt relativistic effective core potential ap- 
proximation, developed in 1970s and applied successfully to many heavy- 
metal compounds (e.g. [10]), to reduce the enormous calculation amount. 
Energy-consistent RECPs dehvered by Stuttgart /Dresden group [16] are uti- 
lized, where 46 core electrons of La are replaced by effective potentials de- 
rived from quasi- relativistic atomic wave functions. Only the outer 11 electrons 
5s^5p^5d^6s^ need to be treated exphcitly in the calculation. Valence basis sets 
are contracted as (8s7p6d4f2g)/[6s5p4d4f2g] including diffuse functions Islpld 
and polarization functions 4f2g[16]. The exponents of the diffuse functions are 
determined according to even-tempered principle, and 4f2g's exponents are 
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4.5, 1.5, 0.5, 0.167, 1.5 and 0.5 respectively. For CI atom Dunning's corre- 
lation consistent triple-zeta basis sets augmented with diffuse functions, i.e. 
aug-cc-pVTZ, are adopted. We use hybrid density functional B3LYP, MoUer- 
Plesset correlation energy correction MP2 and quadratic configuration inter- 
action method QCISD(T) to account for electron correlation effects. Basis 
Set Superposition Error (ESSE) is corrected by Boys and Bernardi's Counter- 
Poise procedure[17]. All the calculations are completed on SGI 1450 server 
with 4GB memory by Gaussian 98 package[18]. 

3 Equilibrium Geometry 

The possible electronic state for the ground state of LaCl may be the singlet 
^T,'^ or the triplet ^A. These two states are both optimized using B3LYP of 
Density-Functional Theory (DFT). The bond length is 0.2520 nm and 0.2576 
nm for the singlet and triplet respectively, and the corresponding experimental 
values obtained from spectroscopic constant are 0.2498 nm and 0.2547 
nm. With optimized geometries QCISD(T) is used to perform high accuracy 
energy calculation. The calculated energy of singlet is about 0.15 eV lower 
than that of triplet, so the ground electronic state is determined to be the 
singlet. This is in accordance with the assignment of [8] although there was 
no compellent experimental evidence in that paper. Theoretical calculations 
here predict LaCl has analogic ground electronic configuration with LaH and 
LaF, which is still to be confirmed by future experiments. 

Table 1 are bond length, bond dissociation energy and harmonic frequency 
obtained with DFT and post-HF methods. Compared with experimental val- 
ues, the bond length is 2.2 pm, 1.6 pm larger at B3LYP, QCISD(T) levels 
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and 0.6 pm smaller at MP2 level. The bias for harmonic frequency is about 
10 cm~^, 1 cm"^ and 12 cm~^ respectively. With respect to bond dissociation 
energy, the values of B3LYP and QCISD(T) nearly equal (4.982 eV and 5.045 
eV), but MP2 value is about 1.8 eV larger. Though by far there is no experi- 
mental report on the dissociation energy, our calculations estimate it is near 
5.0 eV. This estimation is a little smaller than 5.41 eV, predicted by Kaledin 
et al. [19] with a statistic argument that thermochemical properties of the 
lanthanide monohalides may have simple linear relationships. Since B3LYP 
and QCISD(T) can provide more perfect results of energy, we employ them 
to calculate the potential energy curve for the ground state of LaCl. 

4 Potential Energy Curve and Analytical Function 

In order to calculate potential energy curve, first all the possible electronic 
states and dissociation limits of LaCl should be determined based on atomic 
and molecular reaction statics[20,21]. Firstly the separated atomic group me- 
thod can be used to derive the possible electronic state for the ground state. 
The ground electronic states of La and CI arc "^Dg and respectively, both 
belong to SU(n) group. The symmetry is degraded to Coov when La and CI 
approach and form the molecule LaCl. So representations '^Dg and of 
group SU(n) can be resolved into direct sum of the representations of group 

'Dg = e 'Ug e 'Ag, ^ p„ = ^ © ^n.. 

The direct product of atomic representations '^Dg and is reduced as 

(2E+©%©2^J®(2s+©2n„) = ^'=^E+(2)®^'3E-©^'3n(3)©^'^A(2)©^'^$. 
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Obviously the possible electronic states for LaCl are ^'^T,'^, ^'^T,~, ^'^U, ^'^A, 
Secondly the calculated electronic configuration of LaCl is cr^cr^u^Tr^cr^a^ 
TT^cr^cr^TT^cr^, which belongs to the state ^S"*". Therefore the ground state of 
LaCl is X^E+. According to the principle of microscopic reversibihty, the 

dissociation limit for the ground state can be expressed as 

LaCl — ^ La + CI 

Potential energy scans are carried out every 5 pm for a calculation point 
while the internuclear distance is between 0.18 nm and 0.8 nm. Potential 
energy curves we got are shown in Fig 1. Murrell-Sorbie (M-S) function[21] is 
one of the best analytical functions which can describe both short-range and 
long-range interaction very well. Its form reads 

V{R) = -D{1 + aiC + asf + a3C^)exp(-aiO + D, 

where ^ = R/ Re — I, R and Rg are internuclear distance and its equilibrium 
value, D is bond dissociation energy. In addition M-S function can be expanded 
to the following approximate power series at equilibrium position as what 
Dunham did in 1932 [22], i.e. 

The coefficients ai, • ■ • , oe shown in table 2, are calculated by a nonlinear least- 
square fitting of the data of potential energy curve. The formulae to obtain 
force constants from M-S analytical function are /2 = D{a\ — 202), /a — 
2af -|-6aia2 — 603) and = D{?>a\ — 12a\a2 + 2Aaia^). Standard Dunham 
treatment [22] is adopted to calculate vibrotational spectroscopic constants, 
which are listed in table 3 accompanied with those of rotational analyses from 
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experimental spectra. 

ESSE is about 0.12 eV for QCISD(T) method with large basis sets in this 
paper. Considering ESSE correction is often necessary in the high accuracy 
calculation of global potential energy curve. In Fig 1 potential energy curves 
obtained by E3LYP and QCISD(T) are almost the same in the range from 
0.18 nm to 0.4 nm. From 0.4 nm two curves begin to separate. The difference 
is about 0.5 eV when the internuclear distance is 0.8 nm. Comparing with 
the value of bond dissociation energy in table 1, we can find potential energy 
curve got by QCISD(T) does not meet the dissociation limit of ground state, 
although ESSE correction has been made. This problem has been referred to 
by Wang et al.[23]. They argued that the denominator of perturbation energy 
term of single- reference ab initio methods such as QCISD(T) and CCSD(T) 
is connected with the energy gap between Highest Occupied Molecuar Or- 
bit (HOMO) and Lowest Unoccupied Molecular Orbit (LUMO). It was found 
that HOMO-LUMO gap becomes small when the bond separation is far from 
equilibrium, so that the perturbation energy is somewhat large and finally re- 
sults in the too high potential energy. Another possible reason is that RECPs 
and the corresponding valence basis sets are optimized at equilibrium geome- 
try. Multi-reference configuration interaction or multi-reference perturbation 
theory is expected to be efficiently applied to such system far from Rg. 

All the parameters for M-S analytical function and Dunham expansion are 
given in table 2. The fitted dissociation energy is 4.956 eV for E3LYP method, 
which is very close to the value 4.982 eV calculated at equifibrium position. 
Since potential energy of large distance is 0.5 eV above dissociation limit, 
we choose the data of QCISD(T) oi R < 0.4 nm to fit the M-S analytical 
function. Thus the fitted dissociation energy 5.105 eV becomes very close to 
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5.045 eV in table 1. Using the analytical function, most of the calculated 
vibrotational constants except oJeUe, Pe are in rather good agreement with 
those measured in thermal emission spectroscopy experiments [8]. For example 
with B3LYP method, the vibrational frequency LOe for the ground state X^E+ 
is 338.7 cm^^, and experimental value is 341.6 cm~^. The discrepancy with 
harmonic frequency 332.0 cm~^ of table 1 is because of disharmonic effect of 
electron vibrations. Two other spectroscopic constants cUgZe and that the 
experiment didn't analyze out can be also given on theory, the values of B3LYP 
coincide well with those of QCISD(T). With these spectroscopic constants, 
the vibrotational energy levels for the ground state of molecule LaCl may be 
expressed as [cUg — uJeXei^ + + 0.5) + [B^ — aj^v + 0.5)] J( J + 1) — [-De + 

/3e(i^ + 0.5)]J^(J + l)^ + ife-'^^(-'^ + l)^, where v and J correspond to vibration 
and rotation quantum number respectively. 

5 Summciry and Discussion 

The electronic state and dissociation limit of the ground state of lanthanum 
monochloride have been deduced according to the principles of atomic and 
molecular reaction statics and quantum chemical calculations. Molecular prop- 
erties at equilibrium position have also been firstly calculated on B3LYP, MP2 
and QCISD(T) levels with relativistic effective core potential approximation. 
Our theoretical calculations show LaCl has analogic ground state X^E+ with 
LaH and LaF, which is to be confirmed by future experiments. Potential en- 
ergy curves and their M-S analytical functions and Dunham expansions have 
been obtained with B3LYP and QCISD(T) methods. The calculated vibro- 
tational spectroscopic constants from analytical functions are in very good 
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agreement with experiment results. 

All the quantum chemical calculations demonstrate DFT and QCISD(T) 
under energy-consistent RECPs approximation are efficient and reliable to 
evaluate relativistic and electron correlation effects. If consider only bond dis- 
sociation energy, the results of B3LYP and QCISD(T) are totally better than 
those of MP2. We estimate the dissociation energy of gound state is about 
5.0 eV, smaller than the prediction of [19]. Accurate experimental value is 
expected. The potential energy curve of QCISD(T) at large diatomic distance 
is 0.5 eV above dissociation energy, so it doesn't satisfy the requirement of 
dissociation limit. We propose the scheme that fit the data of DFT or the data 
around equihbrium position of QCISD(T) to derive M-S analytical function 
for potential energy curve. According to this scheme, the calculated analytical 
functions of B3LYP and QCISD(T) are both consistent with the dissociation 
limit of the ground state of LaCl molecule. 

The agreement of spectroscopic constants of theoretical calculations and 
experimental measurements shows that Murrell-Sorbie potential function we 
got is a good analytical function to express the potential energy curve of 
lanthanum monochloride. The analytical function can predict transitional fre- 
quency and the intensity in the fine structure of experimental spectra. Its 
application to the investigation of detail collision reaction dynamic processes 
is expected when LaCl is used as catalyst. Thus the analytical function is 
of much realistic importance. The calculation scheme that obtain potential 
energy curve under RECPs approximation and how to use Murrell-Sorbie an- 
alytical function to fit the potential energy curve can be extended to other 
transition heavy-metal compounds. 
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Table 1 

Geometry, bond dissociation energy and harmonic frequency of LaCl 



Method 


Re(nm) 


D(eV) 


t/(cm ^) 


Method 


i?e (nm) 


D(eV) 


1^(0111 ^) 


B3LYP 


0.2520 


4.982 


332.0 


QCISD(T) 


0.2514 


5.045 


342.7 


MP2 


0.2492 


6.788 


353.3 


Exp. 


0.2498^ 




341.6'' 



^ Converted from experimental Be of [8] . 

^ Experimental vibration frequency taken from [8]. 

Table 2 



Parameters and force constants* of analytical potential function for the ground state 
X^S+ of LaCl 

M-S D{eY) ai a2 as /a /s A 

B3LYP 4.956 4.019 0.5252 4.916 0.189 -7.279 227.51 

QCISD(T) 5.105 4.025 0.6610 4.980 0.192 -7.426 233.25 

Dunham ao(eV) ai 02 03 04 05 

B3LYP 37.958 -3.218 6.194 -8.573 9.151 -7.870 5.630 

QCISD(T) 37.970 -3.233 6.384 -9.189 10.271 -9.275 6.974 

* Units for force constants /2, /s, fi are fj-nm~^, fj-nm~^, fj-nm~^ respectively. 

Table 3 

Vibrotational constants* for the ground state of LaCl 

Item We lO^Se t^eXe W^OJeVe lO^Qe 10"/3e 10*De lO^We^e lO^^He 

B3LYP 338.7 9.503 0.9588 7.88 3.579 5.11 2.991 2.72 -2.00 
QCISD(T) 342.0 9.546 0.9571 8.35 3.571 4.76 2.975 2.52 -2.16 

Exp.'" 341.6 9.670 0.9797 7.46 3.640 6.91 3.119 

* Unit is cm~^. 

^ Taken from [8]. 
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X B3LYP point 
* QCISD(T) point 
dash: B.3LYP fitting line 
solid: QCISD(T) fitting line 




Fig. 1. Potential energy curves for the ground state X^S+ of LaCl. The dash and 
sohd Hues are respectively fitting lines of B3LYP and QCISD(T). 
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